Resonances in ID disordered systems: 
localization of energy and resonant transmission 



cn 

O 

o 

(N 



C/3 



o 

(N 
> 

cn 
in 

o 
m 
o 



I 

o 
o 



X 
J3 



K.Yu.Bliokh 

Institute of Radio Astronomy, Kharkov 61002, Ukraine 
Yu.P.Bliokh 

Department of Physics, Technion, Haifa 32000, Israel 
V.Freilikher 

Department of Physics, Bar-Ilan University, Ramat-Gan 52900, Israel 

Localized states in one-dimensional open disordered systems and their connection to the internal 
structure of random samples have been studied. It is shown that the localization of energy and 
anomalously high transmission associated with these states are due to the existence inside the 
sample of a transparent (for a given resonant frequency) segment with the minimal size of order of 
the localization length. A mapping of the stochastic scattering problem in hand onto a deterministic 
quantum problem is developed. It is shown that there is no one-to-one correspondence between the 
localization and high transparency: only small part of localized modes provides the transmission 
coefficient close to one. The maximal transmission is provided by the modes that are localized in 
the center, while the highest energy concentration takes place in cavities shifted towards the input. 
An algorithm is proposed to estimate the position of an effective resonant cavity and its pumping 
rate by measuring the resonant transmission coefficient. The validity of the analytical results have 
been checked by extensive numerical simulations and wavelet analysis. 

PACS numbers: 72.15.Rn, 42.25.Dd, 72.10.Fk 



INTRODUCTION 

Localization of waves in random media has been in- 
vestigated intensively during last few decades. One- 
dimensional (ID) strong localization has received the 
most study, both analytically and numerically. Localiza- 
tion of eigen states in closed ID disordered systems and 
exponentially small (with respect to the length) trans- 
parency of open systems with ID disorder have been 
studied with mathematical rigor (see, for example, Ref. 
1 and references therein). The most common physical 
manifestation of these effects is the fact that a thick stack 
of transparencies reflects light as a rather good mirror. 2J 
Much less evident (though also well known theoretically) 
is that for each (sufficiently long ) disordered ID sample, 
there exists a quasi-discrete set of frequencies that travel 
through the sample unreflected, i.e. with the transmis- 
sion coefficient close to one.^~^ Values of these frequen- 
cies are random, determined exclusively by the structure 
of the realization, and the whole set of frequencies con- 
stitutes a fingerprint of a given random sample. 

High transparency is always accompanied by a large 
concentration (localization) of energy in randomly ar- 
ranged points inside the system. Therefore, an open ran- 
dom ID sample can be considered as a resonator with 
a set of modes (resonances) with high quality factors. 
An important feature of such systems is that, in distinc- 
tion to a regular resonator whose modes occupy all inner 
space, in a ID random sample each mode (frequency) 
is localized inside its own effective 'cavity' whose posi- 



tion is random and size is much smaller than the size 
of the sample. The resonances are sharp and well pro- 
nounced in the sense that their widths are much smaller 
than the distances between them, both in the frequency 
domain and in real space. Therefore, it is of little won- 
der that resonances caused by one-dimensional disorder 
have recently attracted considerable interest of physicists 
looking for the possibility of creating random laser.®" 

In this paper, we investigate the localized states 
(modes, or resonances) in one-dimensional open random 
systems. To understand the nature of the resonances 
in disordered ID systems and to describe them quan- 
titatively, the problem is mapped onto a quantum me- 
chanical problem of tunnelling and resonant transmission 
through an effective two-humped regular potential. With 
the formulas derived on the basis of this mapping, we find 
statistical characteristics of the resonances, namely, their 
spectral density, spectral and spatial widths, as well as 
the amplitudes of the field peaks and transmission coef- 
ficient. Based on these results, one can solve an inverse 
problem, namely, to predict (with some probability) for 
each resonance the position of the effective cavity and 
the maximal pumping amplitude, using the total length 
of the sample and localization length as the fitting pa- 
rameters, and the transmission coefficient as the only 
(measurable) input datum. The analytical results de- 
duced from the quantum-mechanical analogy are in close 
agreement with the results of numerical experiments. 

We study the frequency dependence of the transmis- 
sion coefficient and the spatial intensity distribution in- 
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duced by incident monochromatic waves with different 
frequencies. If the frequency of the incident wave coin- 
cides with one of the 'eigen frequencies' of the sample, 
the energy is genuinely localized in some randomly lo- 
cated area of the size of the localization length. Thus, 
the transparency of this area is abnormally high, whereas 
the transparency of the adjacent segments is exponen- 
tially small. The transparent segment surrounded by 
non-transparent parts of the sample constitutes an ef- 
fective resonant cavity with a high Q-factor. The trans- 
mission coefficient at a resonance is independent of the 
total length, L, of the sample and is determined only 
by the location of the cavity, that is, totally different 
from the typical transmission that decays exponentially 
with increasing L . The maximal transmission is pro- 
vided by the modes that are located in the center, while 
the highest energy concentration takes place in cavities 
shifted towards the input. The number of the resonantly 
transparent frequencies is shown to be L/lioc ^ 1 times 
smaller than the total number of resonances (lioc is the 
localization length). Wavelet analysis of random realiza- 
tions supports the validity of the mapping and provides 
physical insight into the origin of the effective cavities in 
disordered samples. 

FORMULATION OF THE PROBLEM 

We consider a plane monochromatic wave with wave- 
length A and unite amplitude incident from the left 
(x < 0) on a ID disordered sampk^ with randomly fluctu- 
ating refractive index. The transmission coefficient, T(A) 
and the random field amplitude, A = \E{x)\, have been 
studied both analytically and numerically. Typical fre- 
quency dependence of the transmission coefficient of a 
ID randomly layered sample is shown in Fig.l. 

One can sec that along with a continuum of wave- 
lengths for which the transmission coefficient is exponen- 
tially small exp{—2L/lioc)), there is a discrete set of 
points at A-axis (A5, Ac in Fig.l) where T(A) has well pro- 
nounced narrow maxima. The amplitude distributions 
induced by the corresponding waves inside the sample 
are shown in Fig. 2. While at typical frequencies (real- 
izations) the amplitude of the field decrees exponentially 
from the input (Fig. 2, Aa), the resonances, (Fig. 2, At, 
Ac), exhibit essentially non- monotonous spatial distribu- 
tion of the amplitude. 

Important to note that the localization of energy takes 
place for all resonances, not only for those with the trans- 
mission coefficient close to one (Fig. 2, Ab, Ac). The am- 
plitude of a maximum depends on its location in space, 
which in its turn, is uniquely determined by the internal 
structure of the realization. As it is shown below, this 
fact provides a means for evaluating the resonant ampli- 
tude and the coordinate of the point where the resonant 
mode is localized, if the total transmission coefficient is 
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FIG. 1: Transmission coefficient as a function of the wave- 
length. 
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FIG. 2: Amplitude of the field inside the sample as a function 
of the coordinate for three wavelengths marked in Fig. 1. 

known. 

Fig. 3 demonstrates the connection between the ampli- 
tude distribution and the transparency of different parts 
of the sample. The upper picture shows the intensities 
of the field generated by a resonantly transmitted wave 
with T ~ 1 (A5 in Fig. 1) inside the whole sample (black 
curve), and in their left (blue curve), middle (red curve), 
and right (green curve) parts when they are taken as a 
separate sample each. The same dependencies for res- 
onant wave with T < 1 (Ac in Fig. 1) are depicted in 
the lower picture. It is seen that the middle sections 
where the energy is concentrated are almost transpar- 
ent for the wave, while the side parts are practically 
opaque for the resonant frequencies. Wavelet analysis of 
the power spectra of the transparent and non-transparent 
parts presented in Sec. 6 provides physical insight into 
their nature. 

Structures of this type have been studied in the quan- 
tum mechanical problem of the tunnelling and resonant 
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unity) at discrete points corresponding to the eigen lev- 
els of each system. The energy at resonant frequencies 
is localized in a transparent part, and the total trans- 
mission depends drastically on the position of this part. 
More than that, even qualitatively the intensity distribu- 
tions presented in Fig. 3 and the corresponding values 
of the transmission coefficients compare favorably with 
those calculated for a potential profile depicted in Fig. 4 
if the parameters of the effective profile have been chosen 
properly. 



Potential 




FIG. 3: Amplitude of the field as a function of the coordi- 
nate inside the whole sample (black curves), in the left (blue 
curves), middle (red curves), and right (green curves) parts 
taken as a separate sample each. X — Xb and X — Xc. 

transmission of particles trough a regular (non random) 
potential profile consisted of a well bounded at both sides 
by two potential barriers (Fig. 4) 11]. Although the 
physics of the propagation in each system is totally differ- 
ent (interference of the multiply scattered random fields 
in a randomly layered medium, and tunnelling through 
a regular two- humped potential), the similarity of this 
two problems turns out to be rather close. Indeed, 
in both cases the transmission coefficients are exponen- 
tially small for most of frequencies (energies), and have 
well-pronounced resonant maxima (sometimes of order of 



FIG. 4: Two-humped potential profile. 



QUANTUM MECHANICAL FORMULAS 

To map the stochastic classical wave scattering prob- 
lem in hand onto a deterministic quantum mechanical 
one, that is to say, to construct an appropriate effective 
regular potential that provides the same transmission and 
the interior intensity distribution, let us consider an aux- 
iliary problem of the transmission of a quantum parti- 
cle through a potential profile (resonator) consisted of 
two barriers separated by a potential well (Fig. 4). The 
transmission coefficient through the potential can be cal- 
culated in WKB approximation 11], which yields 
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Here Qi — ' ^1 with Ti being the tunneUing trans- 
mission coefficients of the barriers, J/2 is the phase ac- 
quired in the well. 



k{x)dx ^ 2klo 



(2) 



k{x) is the momentum (h = 1), and Iq is the length of the 
well. At 61 = 62 = 6 (symmetrical potential) equation 
(1) transforms into 
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Since 6f 1, it follows from Eqs. (1), (3) that typically 
T <C 1. Resonance transmission (T ~ 1) takes place only 
at the discrete points that correspond to the eigen levels 
of the well: 



J = TT { n 



n = 0,1,2, 



(4) 



Thus, the characteristic spectral interval between the res- 
onances is determined by the requirement A J = tt, or 
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To estimate the half-width of the resonances, 6k, we 
note that in the vicinity of fc^es the transmission coefB- 
cient Eq. (3) takes the form 
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1 + 4^2 (^k-kresYQ^ 

(with Tres = l)j from which it follows that 
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In general (asymmetrical) case Eq. (1) the transmis- 
sion coefficient at a resonance is equal 



T = 

J- res — 
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It is easy to sec that asymmetry of the potential profile 
reduces drastically the resonance transmission: T^es ^ 1 
if 92 <c 92 or 92 ^ 92- In the same approximation the 
amplitude of the wave function in the well is given by 
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APPLICATION TO THE RESONANCES IN 
DISORDERED SYSTEMS 



is, the less is the probability of its occurrence. On the 
other hand, the typical scale in the localized regime is 
lioc- Hence, the minimal and, therefore, the most proba- 
ble size of the effective well is Iq ^ lioc that wc assume to 
be the same in all realizations. Under this assumption, 
different values of the resonant transmission coefficients 
and different intensity distributions at different resonant 
realizations can be reproduced by variations of the loca- 
tion of the well in the corresponding quantum formulas 
presented in Sec. 3 with 



^0 = ho 



(10) 



If the center of the transparent segment of a resonant 

realization is shifted on the distance d from the center of 
the sample, the lengths of the non-transparent parts of 
resonant realizations are 



L - ho 



±d , 



(11) 



with transmission coefficients Ti = exp (— 2/i//;oc)- 
Therefore functions 9, in Eqs. (l)-(9) should be taken 
in the form 



9i,2 = exp 



exp ^ ~ ± — 
^hoc hoc 



(12) 



The above-introduced mapping, Eqs. (10)-(12), al- 
lows the use of Eqs.(l)-(9) for evaluating the correspond- 
ing quantities related to random ID systems. It be- 
comes apparent, for example, that not all eigen modes of 
along random sample provide high transmission (T ~ 1) 
through the system. Indeed, while the distance (in A;- 
space) between eigen modes is 



(13) 



the typical interval between the resonantly transparent 
frequencies obtained from Eq. (5) by substituting Eq. 
(10) is of order 



To use the above-derived formulas for the qualitative 

description of the wave transmission through a random 
sample we have to express the parameters of the ef- 
fective potential through the statistical parameters of 
the disordered system in hand. For an ensemble of 
ID random realizations those parameters are the total 
length L and the localization length hoc that we define 
as hoc = 2L/ < lnT~i Note that in the local- 

ized regime L ^ hoc ^ A. This inequality justifies the 
validity of WKB approximation. To determine properly 
the length of the (effective) well in disordered systems 
we note that the appearance of a transparent segment 
(effective well) inside a random sample is the result of 
a very specific (and therefore low-probable) combination 
of phase relations. Obviously, the longer such segment 
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It is easy to see that Ak/Akres ~ hoc/L <C 1, which 
means that the resonances with T 1 occur much less 
frequently than all other resonances, T < 1. It is phys- 
ically clear, and follows from the fact that a resonance 
with high transmission occurs when the transparent seg- 
ment is located in the middle of the sample (with ac- 
curacy of hoc)^ while the resonances with smaller trans- 
mission are observed in any non-symmetrical potentials 
with arbitrary situated transparent segment. Assuming 
that the transparent segment can be found at any point 
of the sample with equal probability, we infer that res- 
onances with T ^ 1 are encountered L/hoc times more 
rarely than all other ones. 
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It follows from Eq. (7) that typical width of resonances 
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(15) 



that is to say, it decreases exponentially with the length 
of the sample. 

In the same manner, the resonant transmission coeffi- 
cient at resonant frequencies and peak amplitudes can be 
estimated from Eqs. (8) and (9), which after substituting 
Eq. (12) yield: 



e-Kp{2d/lioc) + exp{-2d/lioc) 



(16) 
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8exp{L/lioc - 1 - 2d/lioc) 
[exp{2d/lioc) + cxp{-2d/l. 



(17) 
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Eq. (16) shows that the resonant transmission coeffi- 
cients do not depend on the total length of the sample 
and are governed only by the positions of the areas of lo- 
calization, in contrast to the typical transmission which 
decays exponentially with L increasing. 
If Tres ^ 1, then d~0, and Eq. (17) gives 



1^1 ~ \/2exp {L/2hoc) > 1 



(18) 



When the transparent segment of a random realization 
is shifted from the center, it follows from Eq. (17) that 
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iexp(;i/Z/oc) > 1, if exp(/i///oc) < exp{l2/lioc) 

) {{2I2 - h)/lloc) > 1, if exp(2/2/Z/oc) > exp(/i//;oc)'l > GXp{l2/lloc) 

) {{2I2 - li)/lioc) < 1, if exp{li/lioc) > exp(2/2//ioc) 

I 



(19) 



with li,2 defined in (11). 

It can be shown from Eq. (17) that the maximum of 
the amplitude is reached when the transparent part is 
shifted from the center of the sample towards the input 
on the distance 



d^ --In^ lir,r c^ ~0.27lioc ■ (20) 
This shift is independent of the length of the sample. 



NUMERICAL SIMULATIONS 

To test the validity of the above-introduced analyti- 
cal results, Eqs. (13)-(20), numerical calculations of the 
spectrum of resonances and of the spatial intensity dis- 
tributions at resonant frequencies have been performed 
for more than 10''^ resonances. In the calculations we con- 
sider samples with up to iV = 1000 of layers, and assume 
that the refractive indexes and sizes of the layers are in- 
dependent random variables uniformly distributed in the 
ranges n = 1 ± 0.25 and d — 0.15 ± 0.05 /im respec- 
tively. The wavelength varies in the interval O.bfim < 
X < l.Bfim. 

In accordance with Eq. (14) the average spacing be- 
tween resonantly transparent modes with T ~ 1 depicted 
in Fig. 5 (curve b) is determined by the localization 
length and practically independent on the size L of the 
system. For comparison, the numerically calculated L- 
dependence of the eigen mode spacing is also shown in 
Fig. 5 (curve a). It fits well Eq. (13). Numerically cal- 
culated average half-width of the resonances (Fig. 6) is 
consistent with the analytical expression Eq. (15). The 



localization length had been computed through the trans- 
mission coefficient as koc = —2L < InT > . Note that, 
being a self-averaging quantity Q , koc slightly fluctuates 
from sample to sample, and can be estimated from the 
transmission coefficient at a typical realization. 




FIG. 5: Inverse spacing between eigen modes, (a), and spacing 
between resonantly transparent eigen modes, T ~ 1, (b), as 
functions of the length of the sample. 

Depicted in Fig. 7 is the probability that a mode 
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FIG. 6: Logarithm of the half-width of the resonances as a 
function of the length of the sample. 

localized at a given point provides transmission coeffi- 
cient T (dimensionless coordinate x/L is used). As is 
seen from the picture, the probability of high transitiv- 
ity (T ~ 1) has maximum for modes that are localized 
in the center of the random sample. As the location of 
an effective cavity shifts from the center, the correspond- 
ing transmission coefficient decreases, in agreement with 
Eq. (16). To make the comparison with the analytical 
results more convenient, the same (as in Fig. 7) proba- 
bility density distribution is presented in Fig. 8 as a two- 
dimensional picture where different colors correspond to 
different probabilities. Black line displays the transmis- 
sion coefficient calculated by Eq. (16) as a function of 
the coordinate of the corresponding point of localization. 
Note that shifted to the exit peaks of the field are not 
shown in Fig. 8. They are much weaker, because they 
happen on the background of an exponentially small sig- 
nal. 



FIG. 8: The same probability as in Fig. 7 presented by colors. 
Black line displays the transmission coefficient as a function 
of the coordinate of the point of localization calculated by Eq. 
(16). 

Numerically calculated average normalized field ampli- 
tude of localized modes (pumping rate of effective reso- 
nant cavities) as a function of the coordinate of the cavity 
exponentially increases with the length of the sample in a 
good agreement with Eq. (17). Effective cavities which 
provide the highest pumping rate are shifted from the 
center towards the input in accordance with Eq. (20). 
Interestingly enough, the seemingly rough analogy based 
on only one fitting parameter (localization length) per- 
forms surprisingly well. Indeed, not only the shift is in- 
dependent on the total length and proportional to koc 
as predicted by Eq. (20), but also the coefficient in Eq. 
(20) coincides with that obtained in numerical simula- 
tions with the accuracy 10%. We have also verified that 
the maximal peak amplitudes exponentially increase with 
the length of a sample and by the order of magnitude co- 
incide with the values given by Eqs. (18). 
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FIG. 7: Probability for a mode localized at a given point 
to provide the value of the transmission coefficient T , as a 
function of the dimensionless coordinate x/L. 



WAVELET ANALYSIS OF RANDOM SAMPLES 

Here we discuss briefly why the effective resonant cav- 
ities exist in ID disordered samples, in other words, why 
different segments of a random sample are transparent 
or non-transparent for the wave with given wave num- 
ber k = 2tt/X. It is well known that in the case of a 
weak scattering the so called resonance reflection takes 
place, that is, the reflection coefficient (and, therefore, 
the transparency) of some segment of a random medium 
is determined by the amplitude of the resonance har- 
monic with Qres 

~ 2k in the spatial Fourier transfor- 
mation of its refractive index [l^]. In other words, any 
long enough randomly layered segment, xi < x < X2, 
resonantly reflects as a periodically modulated medium 
whose modulation wave number q is determined by the 
resonance condition q — qres = 2A:, and the amplitude. 
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n^iqres)^ I ^ r?{x)e-'^''''dx . (21) 

J Xx 

This expression is known as window Fourier transform 
with the rectangular window function, and characterizes 
the local spectrum of the investigated function. The rect- 
angular window function has some disadvantages that 
disappear when instead of the window Fourier transform, 
a wavelet transform is applied for the determination of 
the local spatial spectrum of the sample. To investigate 
the local spectra of different parts of the sample we used 
the so-called complex- valued Morlet wavelet We 
first have found the wavelengths for which the trans- 
parency r(A) of the whole sample exceeds 0.5. For these 
wavelengths the spatial distributions of the wave ampli- 
tude and of the amplitude of the corresponding wavelet 
transformation have been compared. The examples are 
shown in Fig. 9. One can see that in the areas where 
field is localized, i. e. in those that should be transparent 
for the given frequency, the wavelet amplitude is strongly 
suppressed, while in the non-transparent parts it has well 
pronounced maxima. Therefore, in accordance with the 
resonant scattering mechanism, the resonant cavity for 
a frequency to arises in that area of a disordered sam- 
ple, where for some reason or other, the harmonic with 
the wave number q^es = 2fc in the power spectrum of 
the (random) refractive index has small amplitude. In 
a sense the wavelet amplitude can be interpreted as an 
effective potential. Note, not all resonance realizations 
exhibit so good correlation with the wavelet analysis as 
shown in Fig. 9. Actually the correlation was observed 
for 70 - 80 % of all resonances. 

CONCLUSIONS 

From the results above it follows that the resonant 
transmission through a disordered ID system occurs due 
to the existence inside the system of a transparent (for 
a given resonant frequency) segment with the size of or- 
der of the localization length. The mode structure of 
the sample, the transmission coefficient, and the inten- 
sity distribution at resonant frequencies depend on the 
positions of the segment. These dependencies are very 
robust, insensitive to the fine structure of the system, and 
therefore can be described by corresponding formulas for 
an effective regular potential profile. The fitting param- 
eters are the total length of the sample and the localiza- 
tion length, which is a self-averaging quantity and can be 
found, for example, from the transmission coefficient at 
typical (non-resonant) realizations. This feature enables 
to estimate the position of an effective resonant cavity by 
measuring the transmission coefficient at a typical and at 
the corresponding resonant frequencies. From the first 
data the localization length can be obtained, the second 
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FIG. 9: Field amplitudes (light lines) and wavelet transforma- 
tions (dark lines) as functions of the coordinate in two random 
realizations. 

one, Tres, could be used to find the asymmetry parame- 
ter d (i.e. the coordinate of the effective cavity) from Eq. 
(16). Then Eq. (17) gives estimation for the pumped in- 
tensity. Therefore, although locally for ID photons there 
is no analog to the quantum mechanical tunnelling (by 
definition the effective energy of ID photons is always 
higher than the effective potential barrier), macroscopi- 
cally (on scales larger than the localization length) the 
problem of the propagation of light through a disordered 
ID system can be reformulated in terms of the effective 
potential profile. 
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